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Abstract 

Measurements of the r lepton polarization and forward-backward polarization asymmetry near 
the TP resonance using the OPAL detector are described. The measurements are based on 
analyses of efg^'r, ^u^Vt, t— > ttVt^ pu^- and r— > aii^^ decays from a sample of 
144,810 e^e~— > T+r" candidates corresponding to an integrated luminosity of 151 pb~^. 
Assuming that the r lepton decays according to V— A theory, we measure the average r po- 
larization near ^/s = Mz to be {Pr) = (—14.10 ± 0.73 it 0.55)% and the r polarization 
forward-backward asymmetry to be A™ = (—10.55 it 0.76 it 0.25)%, where the first error 
is statistical and the second systematic. Taking into account the small effects of the photon 
propagator, photon-Z^ interference and photonic radiative corrections, these results can be 
expressed in terms of the lepton neutral current asymmetry parameters: 

Ar = 0.1456 ±0.0076 ±0.0057, 
Ae = 0.1454 ±0.0108 ±0.0036. 

These measurements are consistent with the hypothesis of lepton universality and combine to 
give Ai = 0.1455 ± 0.0073. Within the context of the Standard Model this combined result 
corresponds to sin^ 9^^^^ = 0.23172 ± 0.00092. Combing these results with those from the other 
OPAL neutral current measurements yields a value of sin^ 9^^^^ = 0.23211 ± 0.00068. 
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1 Introduction 



Parity violation in the weak neutral current results in a polarization of final-state fermion— anti- 
fermion pairs produced in TP decay with the r lepton being the only fundamental fermion whose 
polarization is experimentally accessible using the detectors at the LEP e"'"e~ collider. Consequently, 
measurements of the r polarization provide a means for determining neutral current asymmetry 
parameters which depend on the neutral current vector and axial vector coupling constants. In the 
Standard Model the effective electroweak mixing angle is determined from these couplings, therefore 
the polarization measurements yield a value of sm^ d^^^^ . 

The r polarization, P^, is defined as Pr = (cr+ — cr_)/(cr_|_ + (j_), where cr+(_) represents the cross- 
section for producing positive (negative) helicity lepton^. Assuming unpolarized e+e~ beams, 
the spin-1 nature of the intermediate state implies that the lowest order differential cross-sections 
for cr+ and cr_ can be expressed as: 



(y total dcOS 9r- 

1 

CTtotal dcOS dr- 



^[(1 + (P.) )(1 + cos2^.-) + |(Afb +A™ )cos^.-] 
^[(1 - (P,) )(l + cos2 e^-) + |(ApB - A™ )cose,-] 



where 9^— is the angle between the e beam and the final-state r ; (Ttotai = + ^-]-i<cos6i _<i' 
(Pr) is the average r polarization. 



(Pr) 



i^+J-l<cos6»^_<l i^-J-l<cos6»^_<l 



total 

ApB is the forward-backward asymmetry of the r-pairs, 

ApB = — ; 

0" total 

and A™ is the forward-backward r polarization asymmetry, 

. FB _ [^+ ~ '^-]cose^->0 ~ [^+ ~ ^-]cose^-<0 

O total 



Equation implies a simple dependence of Pr on cos 9^-- 

(2) 



apie aepenaence oi r-^ on cos u^-- '■ 

(Pr) (l + cos^g.-) + |Agg cos^.- 

^^"""^ ""^ (l + COs2^rO + |AFB COs9r- 



where the three parameters, (P^) , A™ and ApB are extracted from data. These parameters 
include contributions from Z° exchange, photon exchange and photon-Z° interference, as well as 
from photonic radiative corrections. At LEP 1, where -^s is near Mz, the pure Z° exchange term 
dominates the polarization. Nonetheless, the small contributions arising from the other components 
must still be taken into account in the interpretation of these parameters in the context of the neutral 
current couplings. The effects of the y/s dependence of the relative sizes of these contributions, 
which are non-negligible near the Z° pole, must also be considered for a precise neutral current 
interpretation. 

"'^By convention, Pr = Pr- and since, to a very good approximation, the and r+ in a given event have opposite 
helicities at LEP: Pr- = -P^+. 



When only the pure Z exchange is considered, the interpretation of the measured values of 
(Pr) and A™ in terms of the neutral current asymmetry parameters At and Ae is remarkably 
simple: 

{Pr) = -Ar and A™ = -^Ae (3) 
where the asymmetry parameters are defined as: 

' " 1 + {9ve/9Aiy ^ ^ 

and the symbols gy^ and qai represent the effective neutral current vector and axial vector couplings 
for lepton i as defined in reference 

The inequality of the Z" coupling to left-handed and right-handed initial-state electrons results 
in a polarization of the Z° which manifests itself as A™ , whereas (Pr) expresses the inequality 
of the Z° coupling to the two chiral states of the final state r leptons. Therefore, once the small 
effects of photon exchange, photon-Z" interference and photonic radiative corrections are taken 
into account, the measurement of (Pr) is directly related to the ratio of the vector to axial vector 
coupling constants for r leptons and that of A™ to the ratio for electrons p|. Consequently, these 
measurements test the hypothesis of lepton universality in the neutral current. 

In the context of the Standard Model, gyi and g^g are related to the effective electroweak mixing 
angle by: 

9ye/9Ae = 1-4 sin^ e^^i\ (5) 

Therefore (Pr) and A™ also provide a precision determination of sin^^^g'*. 

This paper describes a measurement of {Pr) and A™ using the full data sample collected 
with the OPAL detector at LEP during the period 1990-1995 which corresponds to an integrated 
luminosity of 151 pb^^. The OPAL measurement of Aps , which is used as input to this analysis, is 
described in reference|0]. A sample of 144,810 e"'"e~— r+r^ candidate events contained within a 
fiducial acceptance of | cos 9^- \ < 0.90 is used in the analysis. Most of the selected events (88%) were 
recorded with the centre-of-mass energy (y/s) at the peak of the Z" resonance and the remainder, 
referred to as 'off-peak data', were recorded at several distinct a/s values within 3 GeV above and 
below the peak. These new results supersede the measurements reported in Reference 0, which 
were based on an analysis of a 1990-1994 OPAL data sample that was restricted to the central 
region of the detector: | cos6't— | < 0.68. 

The OPAL detector, described in detail in reference ||^, consists of a cylindrical magnetic spec- 
trometer embedded in an electromagnetic (ECAL) calorimeter with presampler and a hadronic 
(HCAL) calorimeter which in turn are enclosed by muon detectors (MUON). A silicon micro- vertex 
detector, vertex chamber, large volume jet chamber, and z-chambers comprise the central tracking 
detector (CT). These are contained in a 0.435 T solenoid with magnetic field aligned along the 
beam axis which defines the z-axis of the detector. The ECAL, consisting of a barrel and two 
endcap arrays of lead glass blocks, provides complete azimuthal coverage within a polar angle range 
of |cos6'| < 0.984. The HCAL and muon detectors are collectively referred to as the 'outer de- 
tectors'. The detector covers nearly the entire solid angle and provides full trigger efficiency!^] for 
e"'"e~^ r+r~ events within the acceptance of this analysis. 

The eue^r, t—>- [iv^Vr, t—>- nUr, T— K z/^, nn^ u^, r— >• K 7r° and t—> 37r^ decays, 
representing a combined branching fraction of 83%, are selected and their kinematic properties used 
to measure the polarization. As nothing in the selection discriminates between charged pions and 
kaons, the r— > ttz/^ and r— * K z/^ modes are analysed as a single channel, denoted as r— > vrz/,-. The 
7r7r° z/^ and r— > K 7r° z/^ modes are dominated by the p^— > vr^vr^ and K*^— K^vr'' resonances, 
respectively, and are similarly treated as a combined channel denoted as pUr- As the 37r^ z/^ 



decays are dominated by the three-prong ai resonance, this channel is labelled as aiu-r- The 
r— > KTCTcUr and r— > K KnUr decay modes are treated as background to the t—>- Stt"^ Ur channel. 
The selection criteria for all channels are based on a likelihood technique and have been optimized 
to reduce the combined statistical and systematic errors. 

In order to perform the polarization measurement, the data are fit with linear combinations of 
left-handed and right-handed Monte Carlo simulated kinematic spectra for each of the five channels. 
The five decay modes do not all have the same sensitivity to the r polarization. For the simplest 
case, the two-body decay of a r lepton to a spin-zero vr meson and r neutrino, r— nUr, the maximum 
sensitivity is provided by the energy spectrum of the vr. The pure V-A charged current decay of the 
r together with angular momentum conservation produces a vr momentum preferentially aligned 
with the helicity of the r. In the lab frame, this means that a vr" produced from the decay of a right- 
handed T~ will, on average, be more energetic than a tt^ produced from the decay of a left-handed 
T~ . For the leptonic channels, r— > ez/gi/^ and r— > fiUf^Ur, the additional unobserved neutrino causes 
a substantial reduction in polarization sensitivity and the decay of a right-handed r~ produces, on 
average, a lower energy charged daughter than a left-handed r~. The r— > pz/,- channel is complicated 
by the fact that the spin-1 p may be produced in two different spin states, the third state being 
forbidden by angular momentum conservation. The longitudinal state dominates the decay width 
and is kinematically equivalent to the vr channel. In addition there is a significant contribution 
from the allowed transverse state which is also sensitive to the polarization but in a manner which 
reduces the polarization sensitivity if only the total visible energy is measured. Much sensitivity 
is regained by spin-analyzing the ti^tt^ decay through measurements of the kinematics of the 
final-state vr" and tt'^I^. The r— aiz/^ channel, which also involves a spin-1 hadron, is similar to 
the p channel except that the ai decays to three pions. As with the r— » pUr channel, sensitivity is 
optimized by spin-analyzing the ai— > tttttt decay via measurements of angles between and momenta 
of the final state pions. The details of the specific variables used for each channel are described 
Section H. 

The maximum sensitivity for each decay mode, defined as l/\fN a where a is the statistical error 
on the polarization measurement using events, is given in Table |l] for Pr=0.0, assuming that all 
the information in the decay, apart from the r direction, is used with full efficiency. A measure of 
the weight with which a given decay mode ideally contributes to the overall measurement of the 
polarization is given by that decay mode's sensitivity squared multiplied by its branching ratio. 
Normalized ideal weights for each decay mode, which are calculated assuming maximum sensitivity 
and perfect identification efficiency, are also given in Table 0. As can be seen, the pUr and 
r— *■ ttz/t- channels are expected to dominate the combined polarization measurement. The actual 
sensitivity achieved in the experiment for the selected event sample is degraded mainly because of 
inefficiencies in the process of selecting a sample of decays and by the presence of background in 
the sample. 

The parameters (P^) and A™ are extracted from the data using a global maximum likelihood 
fit where the data are described by linear combinations of positive and negative helicity distributions 
in observables appropriate to each r decay channel and in 6r-- These distributions are obtained 
from a Monte Carlo simulation in which the r lepton is assumed to decay according to V-A theory. 
For those events in which both r decays have been classified, the analysis explicitly takes into 
account the r+-r^ longitudinal spin correlation by analysing the event as a whole. In so doing, this 
also accounts for experimental correlations between the polarization observables introduced by the 
r-pair selection and decay mode identification criteria. A beneficial aspect of performing a global 
fit to all decay modes simultaneously is the automatic inclusion of the correlations between the 
systematic errors of the different decay modes. 



^The term 'prong' refers to the track of a charged particle originating from the r decay vertex. 
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Branching ratio 


0.18 


0.17 


0.12 


0.25 


0.09 


Maximum sensitivity 


0.22 


0.22 


0.58 


0.49 


0.45 


Normalized ideal weight 


0.06 


0.06 


0.30 


0.44 


0.13 



Table 1: The branching ratios, maximum sensitivity and normalized ideal weight for the five decay 
modes used in the analysis. The ideal weight is calculated as the product of the branching ratio and 
the square of the maximum sensitivity. Presented in the last line of the table is the ideal weight for 
each channel divided by the sum of the ideal weights of the five channels. 



2 Tau Pair Sample 

The selection of the particular decays used for the measurement can be viewed as a two stage 
process where in the first stage a sample of e"'"e~— >■ t^t~ candidates is selected from which samples 
of r— eVeT-'n T— >■ P^i^^i^^t, ttz/^, t— i> pu^- and the three-prong r— aiz/^ decays are identified in the 
second stage. The r-pair selection employs the same basic criteria used in our earlier publications 
for the barrel region of the detector |P, H, ^ with slight modifications to allow for an extension of 
the acceptance into the endcap regions. These are described in detail, along with all aspects of this 
analysis, in Reference |jlO|- The general strategy is to identify events characterized by a pair of back- 



to-back, narrow jets with low particle multiplicity (r-jet). Background from two-photon processes 
is then suppressed by requiring that the events have a minimum total visible energy and significant 
missing transverse momentum when the total energy in the event is low. After removing cosmic ray 
backgrounds, the events which remain are almost entirely lepton-pairs. Events with high measured 
energy that are consistent with being e"'"e~— *■ e"'"e~ or e+e~— >■ p'^p~ are also removed. The polar 
angle of each r-jet with respect to the direction of the e~ beam, 9jet, is determined using charged 
tracks and clusters of deposited energy in the ECAL. Events are selected if the average of | cos 6*^61 1 
for the two r-jets, | cos^jetl, is less than 0.90. The same | cos^^jetl is used as the estimator for the 
magnitude of cos6't— in the analysis. Using this selection, a sample of 144, 810 events is obtained. 

The contributions to the selected events from various physics processes are estimated using a 
number of Monte Carlo data samples which have approximately ten times the number of events 
as the data. The e"'"e~-^ r+r~ signal and e~^e~—>- p^ p~ background are both modelled using the 



KORALZ Monte Carlo generator with the TAUOLA decay package [11| and the e"'"e — >■ e"'"e back- 



ground is estimated using the BHWIDE generator WM. The residual e"'"e — > qq hadronic back- 



ground is simulated using the JETSET Monte Carlo [13] with parameters tuned to fit the global 



event shape distributions of OPAL data [O]. Contributions from non-resonant t-channel two-photon 



processes (e"'"e — >• e"'"e e"'"e , e+e e"'"e p'^ p ) are estimated using the generator described in 
Reference whilst hadronic two-photon processes are modelled using the generator documented 



in Reference [fH]. The potential effects of four-fermion events, including those with hadronic final 



states, are studied using the generator described in Reference [17|. The response of the OPAL 
detector to the generated particles in each case is modelled using a simulation program |TB| based 
on the GEANT |jl9| package. In all cases, the Monte Carlo and real data are treated in an iden- 
tical manner. Estimates of the e"'"e~— * r+r^ selection efficiencies and purities within the different 
fiducial regions of the analysis are obtained using these Monte Carlo samples and are presented in 
Table 0. 



Region 


cos6' 


Efficiency 


Purity 


Barrel 


cos^ < 0.72 


93 % 


98 % 


Overlap 


0.72 < cos^ < 0.81 


75 % 


98 % 


Endcap 


0.81 <cose < 0.90 


79 % 


97 % 



Table 2: r-pair selection efficiencies and purities for the different fiducial regions of the analysis. 



3 Tau Decay Selection 

Starting with this sample of r-pair candidates, each r-jet in an event is classified as one, and only 
one, of r— > euei^r, t— > ^^VhVt-, t—>- TcUr, T— > or r— > Svr^ Ur, or none of these, using a likelihood- 
based identification procedure. Any event having at least one identified r-jet is then used in the 
polarization analysis. From this procedure 129, 902 events contribute to the measurements out of 
the first-stage selection of 144, 810 events. 

Each of the five decay mode selections start with a mode-dependent set of loose cuts which 
define a sample subsequently processed in a binned likelihood selection[^. The likelihood selection 
procedure uses observables which provide discrimination between various decay channels of the r 
lepton and non-r background. The Monte Carlo simulation provides normalized distributions for a 
set of observables, O,, for each of the decay modes. These are subsequently used to calculate for each 
decay channel j, the likelihood, il{Oi), that the measured Oi would be observed. The likelihood 
that decay mode j produces the measured observables in a given r-jet is obtained from the product 
of the likelihoods: C{j) = Hi^^liOi). In order to select decays from mode k, a cut is applied to its 
relative likelihood, L{k) = C{k)/ J2j C{j). From this definition, L{k) lies between and 1 and the 
value of the cut is chosen to maximize the product of purity and efficiency. The requirement that 
decays have large values of L{k) produces a sample with low background normally at the cost of 
efficiency for selecting mode k decays. In this work, each of the five decay mode selections employ 
different observables and therefore exploit a different set of likelihoods, If a r-jet is classified 

in more than one channel after applying the cut to the likelihood, then it is reclassified into the 
channel having the largest relative likelihood. 

Before applying the likelihood selection, fiducial requirements are imposed to remove the small 
fraction of decays having particles entering regions of the detector which are inadequately modelled 
by the Monte Carlo simulation. Furthermore, the discrimination between channels is enhanced by 
dividing the preseleced sample into a set of subsamples according to the detector | cos 9 jet \ region 
(barrel, overlap or endcap as defined in Table 0) and other r-jet characteristics, such as the number 
of electromagnetic clusters in the jet unassociated with tracks. 

Approximately 25% of r leptons decay to Ur and a p meson, which subsequently decays almost 
exclusively to a charged and neutral pion. Consequently, much of the discrimination between the 
r^ pUr and the other decay modes of the r is achieved by identifying photons from a single 7r° 
decay which form an invariant mass with a single charged track consistent with the p mass, mp. 
The ECAL observables use a 'maximum entropy' clustering algorithm[^ which is well suited to 
identifying photons in a r-jet. A cluster in the barrel (overlap, endcap) region of the detector which 
is not associated with a charged track is referred to as a 'neutral cluster' if it has an energy of at 
least 650 MeV (1.25 GeV, 1.0 GeV). When there is only one neutral cluster present in the r-jet then 
it is considered as a vr^ candidate. If there are at least two neutral clusters then the invariant mass 
of the two most energetic clusters forms an observable, mi 2, and the invariant mass of this object 
with the charged track forms the observable referred to as mp. A third invariant mass, m^ef, is 
reconstructed from all neutral clusters and charged tracks in the r-jet and a fourth, mcharged, is 
reconstucted from all charged tracks in the r-jet. Since 9% of r leptons decay via r— 7r27r° Ur, 



and these can appear as a background in the pv^ channel, a fifth invariant mass, mi_prong, is 
reconstructed from the highest momentum track and the neutral clusters in the r-jet. If there are 
more than four neutral clusters, then only the four highest energy clusters are used. These various 
invariant masses are used as likelihood variables. 

The other observables exploited in the likelihood selections are: information from the specific 
energy loss of the charged track as measured in the central jet chamber (dE/dx); the azimuthal 
angle between the highest momentum charged track and the presampler cluster closest to that track 
(0pres); the number of neutral clusters; the ECAL energy associated with the highest momentum 
track in the T-jet (Ea^s); the ratio of Eass to the momentum of the highest momentum track in the 
r-jet (Eqss/p); the ratio of the ECAL energy of the highest energy cluster to the momentum of the 
track with the highest momentum in the r-jet {Fj^ax/v)'-, the ratio of the ECAL energy measured 
in the r-jet from neutral clusters and clusters associated to tracks to the momentum of the track 
having the highest momentum in the jet (Ejet/p); the energy of the neutral clusters not used in 
the calculation of vHp {Eresid)', the energy of the neutral clusters used to calculate mi, 2 (Ei,2); a 
variable describing how well the track from the central detector matches a track segment in the 
muon detectors (CT-MUON); and information from hits measured in HCAL and MUON. 
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Table 3: Observables employed in the likelihood selections used to classify the different decay modes. 
An 'X' indicates that the observable is used in forming the likelihood distribution for the indicated 
decay mode selection. 

The following subsections provide more detail about the likelihood selections for each of the five 
channels. In particular, the subsections specify for each decay mode: the preselection cuts; the 
subsamples upon which likelihood selections are separately applied within the three cos 9 regions of 
the detector; the observables employed in the likelihood selections, which are summarized in Table ^ 
efficiencies within the fiducial acceptance after tau pair selection; the amount of cross-contamination 
from other r decay modes; and the level of background from non-r sources. The observables 
exploited in the likelihood selections are studied and in general are well described by the Monte 
Carlo simulation. Any differences between the data and simulation of these observables are taken 
into account in the assessment of the systematic uncertainties on the polarization measurements. 



Distributions of two observables characterising each decay mode are plotted in Figures |T](a,b)-|(a,b) 
as representative indications of the level of agreement between data and the Monte Carlo simulation. 
The distributions of the likelihoods formed from all observables used in each of the five selections 
are plotted in Figures |l](c)-H(c). Each of these plots comprise a sum of the various likelihood 
distributions used for each decay mode selection and further illustrate the adequacy of the detector 
modelling. The likelihood cut values are optimized separately for each of the various likelihood 
distributions, and not on the combined distributions plotted here. 

3.1 T — ^ ei>ei^T identification 

The preselection of the r— ez/gZ/^ sample requires, within the r-jet, that there be no more than 
three neutral clusters; that no more than three of nine HCAL layers register activity; that no more 
than three muon chambers register activity; and that there be one or two tracks. If there are two 
tracks, it is assumed that the higher momentum track is associated with the electron from the r 
decay and the lower momentum track is ignored. The observables used to form the likelihoods 
include: dE/dx, Eass/p, Ejet/p? CT-MUON, and mp. The likelihoods in the different regions of 
the detector are formed for a subsample of r-jets with the number of neutral clusters equal to 
zero and for a sample with at least one neutral cluster. Distributions of two of the variables used 
in this selection, E^ss/p, and the pull of the measured dE/dx under an electron hypothesis are 
shown in Figure |I]. The final selection produces a sample of 44,083 candidates with an efficiency 
of 92% after tau pair selection within the fiducial region and a background of 4.6%. Most of the 
cross-contamination from r decays arises from ttUt decays (0.9%) and from r-^ pz/^ decays 
(1.3%) along with contributions from a number of other channels (0.8%). The non-r background is 
estimated to contribute approximately 1.6%. 

3.2 T — > fiv^v-r identification 

The likelihood selection for the r— >• [iffj.Pr decays is applied to those r-jets preselected to have a 
single charged track, no more than two neutral clusters and a measured track momentum greater 
than 4% of the beam energy. This sample is divided into subsamples based on information from the 
outer detectors before forming the likelihoods for the different detector regions. The observables 
used in this selection include: dE/dx, Eqs^, Ea^s/p, 0pres, CT-MUON, information from the HCAL 
and muon chambers, and m^. The distributions of the number of HCAL layers hit in the r-jet and 
the number of muon chambers hit for all selected r— > ^v^j.i't decays for both data and Monte Carlo 
simulation of signal and background are displayed in Figure This selection has an efficiency 
of 87% after tau pair selection within the fiducial acceptance and background of 3.3% which is 
heavily dominated by the -kVj. decays with some components arising from the e"'"e~— ^ f^^f^" 
and two-photon processes. This results in the selection of 41,291 decays. 

3.3 T — ^ TTV-T- identification 

The r— s> nUj. decay preselection requires that the r-jet contains one track and that the ratio of 
the momentum of that track to the beam energy be at least 0.02. Within the different detector 
regions this sample is divided into subsamples with zero, one or more than one neutral cluster, and 
subsamples where at least one HCAL layer records hits or where no HCAL layers record hits. For 
these various samples likelihoods are formed from the following observables: dE/dx, Eass/p, Ejet/p? 
(ppres, the number of neutral clusters, CT-MUON, information from the HCAL and muon chambers, 
mi 2, mp, and mjet- Representative distributions of two of the variables important in this analysis 
are seen in Figure]^ for those r-jets selected as r— >• candidates: the pull of the measured dE/dx 
under a pion hypothesis and the number of neutral clusters. A total of 30,440 r— >^ Trur candidates 
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Figure 1: Distributions of (a) Eass/p and (b) pull distribution of measured dE/dx compared to the 
dE/dx expected under an electron hypothesis for selected et'e^'rcandidates. (c) The distribution of 
likelihoods formed from all observables used in the efeZ^r selection. In each figure the points with error 
bars represent the data, the open histogram the ev^.f-j- expectation from Monte Carlo, the hatched 
histogram the cross-contamination from other r decays and the dark shaded histogram the background 
from non-r sources. The distributions include data from all detector regions and subsamples. 
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Figure 2: Distributions of (a) the number of HCAL layers hit and (b) the number of muon chambers hit 

in the r-^ /i^'^^'r selected jets, (c) The distribution of likelihoods formed from all observables used in the 
/vf/iJ^r selection. In each figure the points with error bars represent the data, the open histogram the 
fiv^i^T expectation from Monte Carlo, the hatched histogram the cross-contamination from other r 

decays and the dark shaded histogram the background from non-r sources. The distributions include data 

from all detector regions and subsamples. 




Figure 3: Distributions of (a) pull distribution of measured dE/dx compared to the dE/dx expected under 
a pion hypothesis and (b) the number of neutral clusters, (c) The distribution of likelihoods formed from all 
observables used in the tti^t selection. In each figure the points with error bars represent the data, the 
open histogram the iriyr expectation from Monte Carlo, the hatched histogram the cross-contamination 
from other r decays and the dark shaded histogram the background from non-r sources. The distributions 
include data from all detector regions and subsamples. 



are selected with an efficiency of 75% after tau pair selection within the fiducial region with a 
background level of 26%. Most of the background arises from the r— > pv-r mode (16%) with the 
next largest contributions arising from r— > pv^Vr (5%) and r-^ aiz/i- (2%). The non-r background 
is estimated to contribute approximately 0.2%. 



3.4 T — ^ pUr identification 

The pv^- decay preselection consists solely of the requirement that the r-jet contains one track. 
This sample is divided into the subsamples with zero, one or more than one neutral cluster and 
subsamples where at least one HCAL layer records hits and where no HCAL layers record hits. A 
number of the observables used to create the likelihood selection are also used for the r-^ -kv^ like- 
lihood: dE/dx, the number of neutral clusters, Ejet/p, CT-MUON, information from the HCAL 
and muon chambers, mi^2, nip, and m.jet- Additional observables used to enhance the r— > pVr se- 
lection include: E^esid, Ei_2, and vni-prong- The distributions of two of the significant observables 
used in this selection are displayed in Figure ^ for those r-jets selected as r— >• pVr candidates: m^ 
and the number of neutral clusters. The analysis selects 67,682 pVj. candidates from the r-pair 
sample. Within the polar-angle acceptance, the efficiency is 73%. The background fraction in the 
r— >• pPt- sample is 29% and consists mainly of r— > 7r(K)>27r°z/^ (18%) and r— > vrz/^ (5%) decays. 
The non-r background is estimated to contribute 0.2%. 



3.5 T — > Sl^Vt identification 

For the r— > aiz^T- selection we restrict ourselves to the three-prong mode, which has a branching 
fraction of 9%. It is assumed that all three-pion decays of the r lepton proceed through the 



ai channel p2[| . The r— > aiz/,- r-jet is required to have three charged tracks, none of which is identified 
as a conversion electron. The major contamination for this mode is from Svr > Ivr^z/^ decays. 
Another potential source of background arises from e"'"e~— >■ e"'"e~ events containing a conversion, 
therefore some of the likelihood variables included for the r— >■ a.iu-r selection are designed to suppress 
this contamination. The likelihoods are formed from: dE/dx, Ea^^/p, Ejet/p, Emax/p, ^1,2, CT- 
MUON, information from the HCAL and muon chambers, mi_prong, and m.charged- The distributions 
of Ejet/p and vncharged are plotted in Figure ^ for those r-jets selected as r— > aiz/^- candidates for 
both the data and Monte Carlo simulation of signal and background. The differences between the 
data and Monte Carlo simulation evident in these plots are adequately described by the systematic 
uncertainties assessed for this analysis. In particular, the differences between the data and Monte 
Carlo distributions of vncharged, Figure i|(b), are well described by the systematic uncertainty ascribed 



to the modelling of the r— * aiz/,- decay from the use of different ai models p3|, [2^ . This modelling 
uncertainty is shown by the lightly shaded region of the histogram in Figure |^(b). The number of 
selected r— aii/^ candidates is 22,161. The selection efficiency is 77% and the background is 25%, 
most of which is from r— Svr > lir^Ur and r— > K 7r"'"7r~ Vr decays (20%) and some of which contains 
residual conversions from r— >^ pPr events. The non-tau background in the sample is estimated to 
be 0.2%. 



3.6 Survey of Pair- identification Classes 

The fitting approach adopted in this analysis lends itself to a classification of the data in which 
events are grouped according to the identified decay mode of each of the two r decays in the event. 
For example, if in a given event one r decay is identified as r^ and the other as r^ ttz/,- then 
this event belongs to one pair-identification class, which is denoted as 'vr e'. There are 20 such 
statistically independent classifications including those in which one of the r decays is identified 
and the other is not. An example of the latter would be the pair-identification class denoted as 
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Figure 4: Distributions of (a) the number of neutral clusters for pur candidates and (b) the recon- 
structed p mass, (c) The distribution of likelihoods formed from all obscrvables used in the pi^r selec- 
tion. In each figure the points with error bars represent the data, the open histogram the pi^r expecta- 
tion from Monte Carlo, the hatched histogram the cross-contamination from other r decays and the dark 
shaded histogram the background from non-r sources. The distributions include data from all detector 
regions and subsamples. 
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Figure 5: Distributions of (a) the ratio of the total ECAL energy in the r-jet to the momentum of the 
track with the highest momentum in the r-jet and (b) the reconstructed ai mass for r— > aiVr candidates, 
(c) The distribution of hkehhoods formed from all observables used in the r— > aiVr selection. In each 
figure the points with error bars represent the data, the open histogram the r— > aiVr expectation from 
Monte Carlo, the hatched histogram the cross-contamination from other r decays and the dark shaded 
histogram the background from non-r sources. The lightly shaded region of the histogram in (b) represents 
the T— > aiz^T- decay modelling uncertainty. The distributions include data from all detector regions and 
subsamples. 



'p nid' in which one of the r decays is identified as pv^ and the other is not identified. The 
numbers of data events for each of the 20 possible event pair-identification classes are listed in 
Table ^. Also shown in the table is the expected number of events from Monte Carlo estimates 
from r and non-r sources, scaled by the absolute luminosity. The good agreement between the 
data and the expected numbers of events in the different pair-identification classes evident from 
this table helps validate the overall efficiencies and purities estimated for each class. 
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Table 4: Number of r-pair events in each pair-identification class is presented as the first number 
in each cell. The expected number of events from Monte Carlo estimates using absolute luminosity 
scaling are shown on the second line. The label 'nid' refers to the case where the r decay was not 
identified. 



4 Global Fitting Method 

For the measurements of (P^) and A™ , the distributions described in Equation cannot be 
directly measured as it is not possible to determine the r helicity on an event-by-event basis. 
Instead, distributions of kinematic variables of the r decay products which depend on the r helicity 
are used. These variables, as well as their distributions, depend on the decay mode analysed^]. 

In general, for each r decay channel, i, the decay distribution depends on a set of kinematic 
variables, For positive helicity states the decay distribution can be expressed as + gi{ii) 

whilst for negative helicity states the distribution is — gi{ii)- Consequently, the measured 

decay distribution depends linearly on the weighting of the two helicity states, [p5| : 

1 r/T- 

T^^ = f^{i^)+^rg^{ii). (6) 

For r— >■ eVei^Ti t—>- pVfj.i'T and -kv^- decays, is one dimensional where the relevant kine- 
matic variable ( the charged particle energy scaled by the beam energy. For 
r— >• eUgUr decays, the energy measured in the ECAL associated with the r-jet (xe) is used, whereas 
for r-^ /UZ/^z/t- and r— > TCUr decays, the energy is determined using the momentum of the charged 
particle measured in the central tracking detector. 

For r— pi^T- decays, three kinematic variables enter into the polarization analysis: 6*, the angle 
of the p momentum relative to the r fiight direction in the r rest frame; ip, the angle of the charged 

■^Note that the distributions are the same for the t+ and provided that {Pr) is taken as the t~ hehcity. 



pion relative to the p flight direction in the p rest frame; and nip, the invariant mass of the charged 
particle under a vr"^ hypothesis and the 7r°. This spin-analysis of the p decay recuperates most of 
the sensitivity which would otherwise be lost if only the charged pion momentum were used, as 
discussed in Section |]. These three variables can be converted into a single optimum variable, cj, 
with no polarization sensitivity loss 0. The variable oj is defined by uj=g / f={Ji+ - R_ )/(R+ + 
R_) where R+ and R_ are the population densities of positive and negative helicity r lepton decays, 
respectively, which are functions of and m.p 0. 

The r— > aiz/,- channel is more complicated because the ai decays into three pions. Six observables 
are used in order to improve the sensitivity in the r— >■ aii^T channel[0: the angle {9*) between the 
ai and r momenta in the r rest frame, the angle (■?/') between the perpendicular to the ai decay 
plane and the ai flight direction in the rest frame of the ai, the angle (7) in the ai rest frame 
between the unlike-sign pion momentum in the ai rest frame and the ai flight direction projected 
into the ai decay plane, the Svr-invariant mass, and the two tt+tt" mass combinations present in 
the af ^ 7r^7r"'"7r~ decay. The distribution of the invariant mass of the three charged particles 
assuming them all to be pions, shown in Figure ^d, demonstrates that agreement between the 
data and simulation of this quantity is reasonable within the uncertainties of the modelling of the 
s^iVt- decay. The Monte Carlo distribution depends on the mass and width of the ai as defined 
within the framework of a particular model of r-^ aiz/^ decay I^S] and allowance in the assignment 
of systematic errors is made for r— aiz/,- model dependence. As with the r— pv^, these observables 
are converted into a single optimum variable, cj, which this time depends on the population densities 
of positive and negative helicity r lepton decays which are functions of the six variables mentioned 
above. 

The joint distributions of the r-pair production and decay can be expressed in the improved 
Born approximation as: 

^Oij 51 [(1 + cos^ Or- + |Afb cos^,-) + (7) 

A=±l 

X{{Pr) (1 + C0S2^^,-) + |A™ COS^,-)] X 

I cos^^^- 1) + \Gi{xi, I cos6^-\)][Fj{xj, \ cos6't— |) + XGj{xj, \ cos6't— |)], 

where cjjj is the cross-section to produce an e^e~— r+r^ event in which one r decays via channel i 
and the other via channel j. The first two lines of Equation ^ refer to the production of the r-pairs 
and the third line to the r decays. The summation over A indicates that the summation is over 
positive and negative helicities. The symbol Xi represents the kinematic variable corresponding to 
channel i: xq, x^, xh, up or u^i- The decay distributions for positive-helicity r leptons are given 
by Fi + Gi whereas the decay distributions for negative-helicity r leptons are given by Ft — Gi. 
Fi and Gi represent functions of Xi and | cos 9^— \ after including the effects of the decay mode 
identification procedure, detector response and radiation. The simulation of the underlying detector 
measurements that go into each of the observables used in the analysis, such as track momentum 
and ECAL cluster energies and positions, is checked and corrected if necessary, using various control 
samples as discussed below. Equation |^ includes the correlation between the decay distributions of 
the two r leptons when analysing events in which both r decay channels are identified. 

A binned maximum likelihood fit is performed to simultaneously extract {Pr) and A™ by 
fitting the linear combination of the positive and negative helicity Monte Carlo distributions to the 
data. The values of Xi, xj and cos^^- for each event are calculated and a histogram, binnedQ in 
Xi, Xj and cos 6',— , is then filled for each ^/s. A value for cos^,— of the event is determined from 
I cos 9 jet I and the sign of the charge of the identified r decay. A separate set of histograms exists 
for each combination of decay channel pairs. If only one r decay is identified, then only bins in Xi 

There are ten bins in cos6't-- and ten bins in each of xe, Xjj,, xtt, top, and wai- 



d cos 9t— dxi dxj 



and cos 9^- are filled. The same procedure is performed for the Monte Carlo with a separate set of 
histograms filled for the positive and negative helicity r leptons binned in Xj, Xj and | cos6'^- 1. This 
provides the product [Fi + \Gi\[Fj + XGj] as a function of | cos6't— | in the Monte Carlo, which uses 
the fact that the detector is symmetric in cos 6',— . As a consequence, the forward and backward 
hemispheres use the same Monte Carlo sample. Therefore, the correlations in the Monte Carlo 
samples result in a reduced Monte Carlo statistical error on A™ . 

The Monte Carlo statistics are taken into account in the likelihood fit in the manner described 
in Reference |2^. In order to identify the contribution to the total error arising from the data 
statistical error only, a second fit is performed which does not take into account the Monte Carlo 
statistical errors. The Monte Carlo statistical error is taken to be the quadratic difference between 
the error from the two fits and quoted as part of the systematic error of the polarization results. 

The effects on the measured polarization arising from misidentified r decays are modelled by the 
Monte Carlo simulation. The helicity dependence of the misidentified decays is automatically taken 
into account in the product [Fi + AG »] [Fj + XGj] . Contributions from the small non-r background 
are estimated using Monte Carlo simulations of distributions in the relevant kinematic variables. 
As there is no helicity dependence in this background, these distributions are added to the linear 
combination of the right-handed and left-handed r decay Monte Carlo distributions to form the 
complete reference distributions used in the fit. 

The fit also depends on ApB for which the measured value in the TP — » r+r^ channel at the 
appropriate a/s is used. Separate distributions for the different values of a/s are used in order to 
account for the ApB dependence but a single fit for (P,-) and A™ is performed. Although there 
are potential dependences of the observables in the analysis on the exact value of a/s at which the 
data were collected, the use of beam-energy normalized observables renders the analysis relatively 
insensitive to such effects. However, in order to further reduce any such dependences, the data 
collected with a/s below 90.7 GeV and above 91.7 GeV are analysed using Monte Carlo samples 
generated at fixed centre-of-mass energies where most of the off-peak data were collected. The 
majority of the off-peak data were collected with values of within 0.2 GeV of the values used in 
the Monte Carlo generation. 



5 Polarization Fit Results 

The results of the global fit are: 

{Pr) = (-14.10 ±0.73 ±0.55)% 

A™ = (-10.55 ±0.76 ±0.25)%, (8) 

where the first error is statistical and the second systematic. The correlation between the two 
parameters, including both statistical and systematic correlations, is less than 0.03. Although 
this result uses data collected over a number of different centre-of-mass energies, to a very good 
approximation it can be treated as though it were all collected at a single effective centre-of-mass 
energy of 91.30 GeV. Whereas the systematic error evaluation is discussed in detail in Sections || 
and 1^, this section will outline a number of studies which validate the internal consistency of the 
results. 

The global fit technique has been checked with independent fits to each channel, the results of 
which are presented in Table |^. The weighted average of {Pt) for these fit results differs slightly from 
that obtained from the global fit. This difference is consistent with expected statistical fluctuations 
when comparing results obtained using the global fit which takes into account correlations and 
a weighted average which does not. The values obtained for the five different channels are also 
consistent with each other and the global fit values. A of 4.9 for four degrees of freedom is found 



when comparing the five values of (Pr) to the value from the global fit and 2.1 for four degrees of 
freedom when comparing the values of A™j . 
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44,083 
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3.3% 
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Koi (%) 


-18.7±2.5 
-8.9±2.6 


-16.3±2.7 
-10.6±2.8 


-13.8±1.2 
-11.5±1.3 


-13.3±1.1 
-10.6±1.1 


-11.6±2.8 
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Table 5: The number of decays in the sample, selection efficiency after tau pair selection within 
the fiducial acceptance and background for each decay mode analysed. Results of independent fits 
for the individual decay modes are also presented where the error quoted represents that arising 
from the data statistics only. The measurements from the individual channels are correlated and 
therefore should not be combined in a simple average. 

An indication of the validity of the Monte Carlo simulation of the kinematic variables used 
in the fit, as well as the assumed efficiencies and purities, is provided by the one- dimensional 
distributions of the relevant kinematic variables for the five channels. The distributions combined 
over all cos 6r- bins are shown in Figure ^ for both the data and Monte Carlo. The comparison 
between data and Monte Carlo simulation for these distributions, shown on each plot, suggests 
that the efficiencies and background are adequately simulated in these distributions for all modes. 
Also shown are the Monte Carlo distributions of the variables for positive and negative helicity r 
lepton decays and their sum including non-r background, assuming the value of (Pr) quoted in 
Equation These plots illustrate the sensitivities of the measured kinematic distributions to the 
polarization. 

Another tool for investigating the internal consistency of the analysis is afforded by comparing 
the fitted values of {Pr) and A™ for each event pair-identification class as defined in Section 
These comparisons are shown in Figures |^ and ^ where the results of the fits both in graphical 
and numeric form for each pair-identification class in the global analysis are presented. The 
probability describing the statistical significance of the different values from the global fit valueQ for 
(Pr) and Ap^i indicate internal consistency. The ideogram formed from the sum of the individual 
Gaussians is overlayed on the data points. This illustrates that for both {Pj.) and A™ the spread 
in fit values for the different pair-identification classes is symmetric about the global fit value and 
the peak is consistent with that value. 

High statistics internal consistency is also examined for events which have both r decays classified 
compared to those where only one r decay is classified. These comparisons, also shown in Figures |^ 
and H, indicate strong overall internal consistency of the results. 

As a further check on the validity of the fit, the results of fits for P^- performed independently 
in ten bins of cos 6r- are shown in Figure For the fit in a particular cos 6r- bin, an expression 
analogous to that shown in Equation ^ is used in which ((Pr) (1 + cos^ 9^--) + fA™ cos^,--) in the 
second line of Equation |^ is replaced by (-Pt(1 + cos^ 9^-- + |Afb cos^^^-))- This substitution uses 
Equation |[ Overlaying these points is a curve which represents the expectation value of Pr as 
a function of cos 6'^- from Equation |^ where the values of the (Pr) and A™ from the global fit 
of Equation ^ are used. The results of the ten independent fits are in good agreement with the 
expectations from the global fit: the is 6.8 for eight degrees of freedom when comparing the ten 

'"'These are calculated using data statitistical errors only. 
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Figure 6: Distributions in the kinematic variables used in the fits as discussed in the text for the r— > eve^n 
Ixv^Vr, T— > irvTi pi^r, and r— > aii'r channels where the data, shown by points with error bars, are 
integrated over the whole cos 6^- range. Overlaying these distributions are Monte Carlo distributions for the 
positive (dotted line) and negative (dashed line) helicity r leptons and for their sum including background, 
assuming a value of {Pr) = —14.10% as reported in the text. The hatched histogram represent the Monte 
Carlo expectations of contributions from cross-contamination from other r decays and the dark shaded 
histogram the background from non-r sources. The level of agreement between the data and Monte Carlo 
distributions is quantified by quoting the and the number of degrees of freedom. 
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Figure 7: Internal consistency of the (P,-) results investigated as a function of the number of r decays 
classified in the event and by pair-identification class. The ideogram formed from the sum of the individual 
Gaussians is superimposed on the pair-identification results. The probabilities of the spreads about the 
global fit value are shown for each subsample and show good internal consistency in all cases. The label 
'nid' refers to the case where the r decay is not identified. 
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Figure 8: Internal consistency of the A™ results investigated as a function of the number of r decays 
classified in the event and by pair-identification class. The ideogram formed from the sum of the individual 
Gaussians is superimposed on the pair-identification class results. The probabilities of the spreads about 
the global fit value are shown for each subsample and show good internal consistency in all cases. The 
label 'nid' refers to the case where the r decay is not identified. 
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Figure 9: Tau polarization, (P-j-) , as a function of cos9^-. The data points represent the Pr values 
obtained from a global fit to all channels in each cos 6^— bin. The error bars represent data statistical 
errors only. The curve represents the expectation from the global fit result: (Pr) = —14.10% 
—10.55% as reported in the text. 



and A™ 



values of Pr to the expected value from the global fit where only the data statistics are included in 
calculating the x^- 

Because a non-negligible amount of data is collected off-peak, measurements of (Pr) and 
A™ are also performed separately for data collected at centre-of-mass energies below, on and 
above the Z° resonance peak. The three statistically independent measurements are quoted in 
Table |^. In order to compare the consistency between these measurements, they are all converted 
into measurements of At or using the technique described in Section ^ These asymmetry val- 
ues, also quoted in Table |^, are consistent with each other. Because the conversions from (Pr) to 
At and from A™ to Ae assume the Standard Model centre-of-mass dependence of {Pt) and A™ , 
the agreement between the At and Ac values in Table ^ indicates that the data are consistent with 
the Standard Model expectations for this centre-of-mass dependence. 



(GeV) 


(Pr) Kol 
(%) (%) 




89.5 ± 0.2 
91.25 ±0.05 
93.0 ±0.2 


-15.9 ±3.3 -8.0 ±3.4 
-13.52 ±0.78 -10.52 ±0.81 
-18.1 ±2.9 -11.3 ±3.0 


0.186 ± 0.039 0.125 ± 0.040 
0.1393 ± 0.0081 0.145 ±0.011 
0.176 ± 0.028 0.146 ± 0.040 



Table 6: Global fit values of (P^) and A™ for data collected below, on and above the Z° resonance 
peak. The luminosity weighted values of ^/s are quoted in the first column where the error reflects 
the spread in ^ values of the data combined in each fit. The neutral current asymmetry parameters 
with their statistical errors, based on the data collected at the different centre-of-mass energies, are 
also quoted. 

6 Detector Related Systematic Errors 

Because the Monte Carlo simulation provides the positive and negative helicity reference distribu- 
tions in the fit and is used for the likehhood selections, it is necessary that the detector response 
be accurately modelled and a systematic error attributed to each relevant discrepency between the 
detector simulation and data. The approach taken in this analysis is to determine the accuracy of 
the detector modelling of specific measured quantities, to vary that quantity in the simulation, and 
then to propogate the influence of that variation through a complete analysis of all decay modes. 

High purity control samples of muons with momenta of approximately 45 GeV from e"'"e~— > 
events are used to determine corrections to the simulation of the momentum scale and resolution of 
the central tracking detector. The systematic uncertainties of these corrections yield a transverse 
momentum scale uncertainty of 0.14% in the barrel region of the detector with higher values of 
0.8% and 0.4% in the overlap and endcap regions, respectively. Studies of measurements of the 
mass of the Kg from Kg— > tt+tt" in e+e"— > qq hadronic events provide a calibration point of the 
momentum scale at the lower energies. The tracking resolution corrections are cross-checked at 
lower energies using the transverse momentum distributions in e^c^^ c^e~fi~^fi~ two-photon pro- 
cesses. Modelling uncertainties of the cos^ measurements of the tracks in the tracking chambers 
are evaluated independently of transverse momentum using back-to-back tracks from e"'"e~— > /^"'"/x" 
events and cross-checked with other detectors such as the presampler. Pure samples of electrons 
with energies of approximately 45 GeV from e'^e"— > e"'"e~ events are used to determine corrections 
to the simulation of the energy scale and resolution of the ECAL. These corrections were evolved 
to lower energies using the ratio of the deposited energy to measured momentum for electrons in 
e"''e~^ e"''e~e^e~ two-photon processes and in high purity ez/gZ^r samples. Uncertainties of 0.3% 
on the ECAL energy scale in the barrel, 0.6% in the overlap and 0.4% in the endcap are estimated 
from these studies. The one standard deviation errors on the energy scale and momentum scale are 
used in assessing the systematic errors on (P^) and Ap^j from an analysis using rescaled energy 
and momenta which takes into account the correlations between channels. In a similar manner, 
systematic errors associated with uncertainties in the parameters used to describe the resolutions 
of the ECAL and tracking detector resolution are also assigned. 

High purity muon and r— > samples are used to correct the modelling of the response of the 
HCAL and muon chambers to muons and hadrons. Variation of the magnitude of these corrections 
is used to assess the systematic error on (P^) and Ap^j associated with this modelling. Correct 
modelling of the dE/dx measurement is achieved by studying the response of the tracking detector 
to high purity r— > ez/gi^r and r— > iiP^Ur samples selected without using dE/dx information. The 
corrections applied to the dE/dx simulation are changed in order to assess the sensitivity of (P^) and 
A™i to this modeUing. The effects of uncertainties in the amount of material in the central detector. 



which potentially affects the photon conversion background in the ai channel, were studied and 
found to have a negligible influence on the polarization measurement. 

The uncertainty in the modelling of the lateral spread of the electromagnetic and hadronic 
showers in the ECAL contributes significantly to the overall systematic error. This is particularly 
relevant for the separation of the r— nUj. and r— > pu^- samples. The influence of these uncertainties 
on the polarization measurement is estimated by varying the thresholds in the cluster definitions 
for the simulation and from the cluster position and position resolution uncertainties. The ECAL 
cluster position resolution is also sensitive to the lateral shower spread in the ECAL. A control 
sample of electrons is used to improve the modelling of the ECAL cluster position resolution and 
to assign uncertainties to this modelling by varying the magnitude of the corrections applied to 
the simulation. Further checks of this class of systematic error were performed by studying the 
stability of the results from the likelihood selection when excluding individual observables related 
to showering. 

There is also a potential systematic error on A™ related to charge mis- assignment, which in 
OPAL is negligible. 

The contributions to the uncertainty on {Pt) and A™ from these various sources are shown in 
Table |^. These are obtained by varying the associated detector-level quantities, such as corrections 
to the simulated momentum or HCAL layers hit, and redoing the analysis from the tau pair selection, 
through the calculation of physics observables such as m^ and decay classification, to the global 
fit for (Pt) and A™ . The table includes the systematic errors for each of the five channels in 
addition to the errors for the global analysis. The systematic correlations between channels are 
automatically accounted for in this analysis and fully incorporated into the systematic error quoted 
for the global fit results. 
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Table 7: Tabulation of systematic errors contributing to (P^) and Ap^j when these asymmetries 
are expressed as a percentage, for each of the five decay modes analysed and the global fit. In each 
column the error on (Pr) is given first followed by that on A™ . Systematic error correlations 
between the five channels are fully incorporated into the systematic error on the global result. In 
the second to sixth columns a dash indicates that the listed effect contributes less than 0.05%. 



7 Physics Related Systematic Errors 



Another class of systematic uncertainties relates to our knowledge of r production and decay. In 
this category are the errors on measured branching ratios of the different r decay modes. The 
branching ratios used are obtained from the averages of the measurements in Reference The 
error on (P^) and A™ associated with the uncertainty of each branching ratio is estimated by 
varying the value used in the global analysis by ±1 standard deviation about its average value. 

The uncertainty in the modelling of the r— > aiz/,- decay introduces systematic errors both in 
the r— > ^iUt- channel and in the r— > pu^- channel where the t—>- a,iVj. decays represent a significant 
fraction of the selected decays. Two contributions to the r— > a,iVj. modelling uncertainty are con- 



sidered: one being the uncertainty in the mass and width of the ai as obtained from Reference |^ 
and the other obtained by comparing two independent theoretical treatments of the r-^ aiz/^ de- 
cay 



23, 24, 30 



In addition to the r— a.iVr modelling uncertainty, the modelling of the r— Svr > Itt^Ut decays 
potentially introduces a further uncertainty in the analysis of the r— aiUj. channel. This is studied 
using the Monte Carlo simulation by varying the ujtt^ contribution to the Svrvr^z/^ mode 

and by flipping the helicity of the r in the t—>- 3/t > ln^Vr decays. The latter effect is found to 
contribute negligibly to the systematic error. 

Two smaller sources of error also fall into this general category of systematic error: the error 
associated with the measured value of ApB for e+e~-^ r+r~ events which is obtained from Ref- 
erence 0] and the uncertainty of the simulation of radiation. Both initial and final state radiation 
are sufficiently well modelled in the Monte Carlo and contribute negligibly to the systematic un- 
certainties. The simulation of radiation in the decay of the r is treated in the same manner as is 
described in Reference [§. The uncertainty in the modelling of radiation in all modes is found to 
contribute negligibly to the overall systematic error. 

The contributions arising from the various systematic errors are summarized in Table |^ for each 
of the independent analyses and for the global analysis which takes into account the correlations 
between channels. 

Because the fits also depend on our knowledge of the non-r sources of background, there are 
potential systematic errors arising from the uncertainty in the production cross-sections. Varying 
the cross-sections within their errors in the reference distributions used in the fits makes negli- 
gible changes to the polarization results. To cross-check the contributions of non-r background, 
distributions of acoplanarity, acolinearity and total event transverse momentum, which would have 
regions enhanced in the two-photon, e"'"e~— >• , and e"'"e~— > e"'"e~ events, were studied and 

show no indication of uncontrolled sources of non-tau background. The samples of two photon and 
e"'"e~— >■ e"'"e~ events are further enhanced by examining off-peak data where, again, there is no 
evidence of problems with these sources. Another cross-check is provided by studing the numbers 
of events in the different pair-identification classes, where, once more, there is no indication that 
additional systematic errors associated with non-tau backgrounds need to be quoted beyond those 
included in accounting for the detector response related systematic errors. 



8 Neutral Current Asymmetry Parameter Results 

The polarization measurements quoted in Equation ^ are consistent with our previous measure- 
ments but with a total error that has been reduced by nearly a factor of two. The results are 
also consistent with those published by other LEP collaborations ||31|, ^ ^ . 

As discussed in the introduction the Standard Model gives predictions for [Pt) and A™ in 
terms of the mass and width of the TP , g^i and qm- The connection is made through the 



neutral current asymmetry parameters as defined in Equation |. The measurements of (Pr) and 
A™ are dominated by tlie propagator but also include small contributions from the photon 
propagator, photon-Z° interference and photonic radiative corrections. ZFITTER provides the 
^/s dependent non-Z" propagator contributions to (Pr) and A™ as well as higher-order corrections 
to Equation ^ within the context of the Standard Model. This allows the measured parameters to 
be expressed in terms of At and Ae- 

Ar = 0.1456 ±0.0076 ±0.0057, 
Ae = 0.1454 ±0.0108 ±0.0036. 

Within the context of the Standard Model these results can be interpreted as measurements of: 

gvr/gAr = 0.0732 ±0.0048, 
gve/gAc = 0.0731 ±0.0057, 

where the statistical and systematic errors of Ar and Ac are added in quadrature before calculating 
the errors on gyr/gAr and gve/gAe- The agreement between these two values indicates that the data 
are consistent with the hypothesis of lepton universality. This test of lepton universality can be 
expressed as a measurement of the ratio of the tau vector to axial vector couplings to the electron 
vector to axial vector couplings: 

^^^^ = 1.00 ±0.10. 

gVe/gAe 

If universality is assumed, these data can be averaged and expressed as a single leptonic asymmetry 
parameter: 

Ae = 0.1455 ±0.0073. 

Using Equations ^ and |, this result can be expressed as: 

sin^^^ffP* = 0.23172 ±0.00092. 

This measurement of sin^ ^gg'* is of similar precision to other individual measurements at LEP 
using various techniques and is in agreement with the value of sin^ 6^^^^ obtained from a Standard 
Model fit to all LEP electroweak data, including previous measurements of the r polarization [p4 |. 
It is also consistent with a determination of sin^ from a measurement of Alr by the SLD 



collaboration 35 



9 Combined Lineshape and Asymmetry Results from OPAL 

The tau polarization results can be combined with the measurements from the leptonic partial 
widths and forward-backward asymmetries published by OPAL[^ to provide measurements of 
gye and g^e for electrons, muons and r leptons separately. While Ar and Ae provide measure- 
ments of gye/gAi for the r and electron, the leptonic partial widths of the provide mea- 
surements of QYi"^ + QAe^ for all three lepton flavours. The forward-backward asymmetries de- 
termine {gyegye) / {gxegAi) for electrons, muons and r leptons yielding additional information about 
gwr/gAr and gve/gAe, and provide the only means of measuring gy^/gAfi at LEP. 

The OPAL measurements of the hadronic and leptonic cross-sections and leptonic forward- 
backward asymmetries at the Z^ pole are summarized in terms of nine model-independent Z^ pa- 
rameters as defined in Reference 0: the mass (Mz), width (Pz) and the hadronic pole cross-section 
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Table 8: The first nine parameters are the result of fitting the model-independent Z" parameters 
to the measured cross-sections and asymmetries measured by OPAL The parameters Ae and 
At are the result of the current analysis of the r polarization. 
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Table 9: Error correlation matrix for the 11 parameters entering the fit for the leptonic couplings 
presented in Table 
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Table 10: The leptonic neutral current asymmetry parameters obtained from a fit to the Z pa- 
rameters given in Table In the last column is given the value of the parameter calculated in the 
context of the Standard Model assuming the parameter variations given in the text. 





Aq 


A^ 




Aq 


1.00 


-0.43 


-0.09 


Af^ 


-0.43 


1.00 


0.04 




-0.09 


0.04 


1.00 



Table 11: Error correlation matrix for the measurements of the leptonic neutral current asymmetry 



parameters, which are presented in Table 10 



(cr°) of the Z° resonance; the ratios of the hadronic to leptonic partial widths {R^.R^ and Rr); 
and the leptonic pole forward-backward asymmetries (Ap^jApg and ^^)- The values of these nine 
parameters, together with the polarization asymmetries reported here, are displayed in Table ^ and 
the error correlation matrix in Table ^. 

Information from the tau polarization and forward-backward asymmety measurements can be 
combined to provide measurements of the three leptonic neutral current asymmetry parameters, 
^c, Afj, and At ■ These and their error correlation maxtrix are shown in Tables |10| and |lT[ The 
results for the three lepton species are consistent with each other and agree well with the prediction 
of the Standard ModelQ which is also shown in Table |10[ Assuming lepton unversality. 



Af = 0.1424 ±0.054. 



This can also be expressed in terms of the effective weak mixing angle: 

sin^^^'^fP' = 0.23211 ± 0.00068. 



The eleven model-independent parameters listed in Table || can also be used to determine 
the neutral current vector and axial-vector couplings for each lepton species. The results and the 
error correlation matrix are given in Tables |12| and O and are illustrated in Figure |10|. Some of 



^The Standard Model calculations require the full specification of the fundamental Standard Model parameters. 
The main parameters are the masses of the Z° boson (Afz), the top quark {mt) and the Higgs boson (toh), and the 
strong and electromagnetic coupling constants, as and a. As in Reference [|| the calculation of Standard Model 
predictions use the following values and ranges: Mz = 91.1856±0.0030 GeV, mt = 175±5 GeV, mn = 150tgo° GeV, 
as = 0.119 ± 0.002, and a{Mz^)^^ = 128.886 ± 0.090 . The choice of these parameter values and ranges is discussed 
in Reference 



the couplings have up to approximately 50% correlations between their errors. Evident from these 
results is the universality of the coupling constants which can be quantified in terms of the ratios 
of the couplings: 



^ = 1.0011 ±0.0025, ^ = 1.0021 ±0.0029, — = 1.0009 ± 0.0027, 

^Ae g^^ 



9Ye 



1 12±°-2'^ 

-'-•-'-^='=0.21' 



9yt 



1 Ofi+0-1^ 



9Yt 



0.94 ±°:}^ 



The axial- vector couplings of the different lepton species are found to be the same at the 0.3% level. 
The errors on the inter-species ratios of the vector couplings are much larger because of the smaller 
size of the vector couplings themselves, but within this reduced sensitivity, again no significant 
differences are observed. These errors have been significantly reduced by adding information from 
the tau polarization to the lineshape and forward-backward asymmetry measurements [^. 

Combining the values of the coupling constants from the different lepton species under the 
assumption of lepton universality yields the values 



9Ae 



-0.50089 ± 0.00045, gve = -0.0358 ± 0.0014. 



where the correlation between qm and gy£ is —19%, which is approximately the same as the —21% 
correlation between and sin^ These are in good agreement with the predictions of the 



Standard Model, which are included in Table 12. 
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Table 12: Axial- vector and vector couplings obtained from a fit to the parameter set given in 
Table ^. In the last column we give the values of the couplings calculated in the context of the 
Standard Model assuming the parameter variations given in the text. 



10 Summary 

Measurements of (Pr) and Ap^j have been made using the complete LEP I data sample of OPAL. 
The results are based on a simultaneous analysis of r— >■ eve^ri t— > f^i^f^i^r, ttz/,-, pu^- and 
r— >■ sliUt- decays from a sample of 144, 810 e"'"e~^ t~^t~ candidates collected over nearly the entire 
solid angle of the OPAL detector. Under the assumption that the r lepton decays according to 
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Figure 10: g\i vs gxi as determined from the OPAL measurements of tlie leptonic partial widths of the 
TP , forward-backward asymmetries and tau polarization measurements. The ellipses represent the 68% 
confidence level contours in the gy£-gA£ plane for each lepton species separately (dotted and dashed) and 
for all leptons assuming universality (solid). The central values are displayed at the centre of the ellipses 
as a circle, square, triangle and star for electrons, muons, tau leptons and all leptons under universality, 
respectively. The Standard Model prediction is shown with variations from the top quark mass (170 to 
180 GeV) and Higgs mass (90 to 1000 GeV) indicated. The OPAL tau polarization measurements of At 
and constrain gy£ and g^e to lie in the shaded region at the 68% confidence level. 
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Table 13: Error correlation matrix for the measurements of the axial vector and vector couplings, 
without assuming lepton universality, which are presented in Table |12. 



V— A theory, the average r polarization near ^/s = Mz is measured to be (Pr) = (—14.10 ± 0.73 ± 
0.55)% and the r polarization forward-backward asymmetry to be A™ = (—10.55 ±0.76 ±0.25)%. 
Taking into account the small effects of the photon propagator, photon-Z° interference and photonic 
radiative corrections, these results can be expressed in terms of the lepton neutral current asymmetry 
parameters: 

Ar = 0.1456 ±0.0076 ±0.0057, 
Ae = 0.1454 ±0.0108 ±0.0036. 

These measurements are consistent with the hypothesis of lepton universality and combine to give 
Ai = 0.1455 ± 0.0073. Within the context of the Standard Model this corresponds to sin^ 6*'*^* = 
0.23172 ±0.00092. 

Combining the information from the tau polarization results with the results of the other OPAL 
neutral current measurements yields values for the vector and axial-vector couplings which are the 
same for all lepton species and gives 

gj^^ = -0.50089 ± 0.00045, gyi = -0.0358 ± 0.0014. 

Expressing these results in terms of the electroweak mixing angle gives 

sin^ = 0.23211 ± 0.00068. 



This is consistent with the Standard Model and with the current world average value 1 34 . 
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